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Abstract 

We present a simple, general energy functional for fer- 
romagnetic materials based upon a local spin density 
extension to the Stoner theory of itinerant ferromag- 
netism. The functional reproduces well available ab 
initio results and experimental interfacial energies for 
grain boundaries in iron. The model shows that inter- 
granular cohesion along symmetric tilt boundaries in 
iron is dependent upon strong magnetic structure at 
the interface, illuminates the mechanisms underlying 
this structure, and provides a simple explanation for 
relaxation of the atomic structure at these bound- 
aries. 

Iron and its steel alloys, exhibit two distinct out- 
standing physical properties, high strength and high 
magnetic response. Little is known about the di- 
rect relationship, if any, between the microscopic ori- 
gins of these properties. The exchange interaction, 
preference for alignment of spins to reduce inter- 
electronic Coulomb repulsion through the statisti- 
cal avoidance of Fermions, drives itinerant ferromag- 
netism. Fundamental topological excitations of the 
crystalline lattice govern mechanical response: dis- 
locations mediate plasticity, grain boundaries con- 
trol microstructure and inter-granular cohesion. In 
this letter we present a new model for itinerant fer- 
romagnetic iron and demonstrate that the exchange 
interaction plays a significant stabilizing role in grain 
boundaries, thereby contributing to the strength of 
inter-granular cohesion. 

The fundamental role played by exchange in con- 
trolling the mechanical response of magnetic mate- 
rials is at present poorly understood for lack of an 
appropriate microscopic theory. Studies in iron to 
date are limited to either full-blown ab initio spin- 
dependent electronic structure calculations or sim- 
ple interatomic potentials such as the embedded atom 
method (EAM)|^. While describing all of the cor- 
rect physics, ab initio calculations treat explicitly too 



many degrees of freedom to allow studies of the com- 
plex structures of all but the simplest extended de- 
fects. For this reason, total energy ab initio studies 
of iron to date have been limited to small clusters en- 
compassing the behavior of the simplest (S = 3(111)) 
grain boundary j|, ^ |^. Interatomic potentials, while 
practical for complex systems, deal with too few de- 
grees of freedom to treat itinerant exchange prop- 
erly. Nonetheless, these potentials have been useful 
in studies of grain boundaries dislocations and 
even fracture |Q. Energies calculated for grain bound- 
aries with the EAM|^, however, are exaggerated by 
about a factor of two when compared with experimen- 
tal values]^, beginning to approach even the experi- 
mental surface energy |lO[. Here, we present a simple, 
general model which gives much better results and 
reveals the mechanisms stabilizing the boundaries. 

In response to the above weakness, Krasko[|llj has 
introduced recently an atomistic potential which in- 
cludes a prescription for estimating the exchange en- 
ergy of each atom from its local environment. Here, 
we follow an alternate, microscopic route. In ac- 
cord with the microscopic origin of itinerant ferro- 
magnetism, we consider the extended nature of the 
electrons and determine only the spin-moment cou- 
pling constant, the Stoner parameter Q , as a func- 
tion of the local atomic environment. 

Microscopic Approach — Such intermediate, elec- 
tronic structure based descriptions of iron have been 
developed in previous studies of iron|l^, [llj, but for 
ideal crystalline systems. Our approach to the study 
of defects is to first identify the smallest, physically 
reasonable set of degrees of freedom from these stud- 
ies, and then extend this treatment to include inho- 
mogeneous systems. 

Hasegawa and Pettifor[^ reproduce the experi- 
mental P — T phase diagram of iron by combining 
a tight-binding description of the d— bands with a 
treatment of spin-fluctuation effects. They note that 
as r — > 0, spin fluctuation effects become unim- 



1 



portant and their spin-fluctuation theory reduces to 
the traditional mean-field Stoner theory of itinerant 
ferromagnetismp^ . They identify the temperature 
at which this happens to be Tf w 500 K. Zhong, 
Overney, and Tomanekjl^ take up the fact that 
mean-field Stoner theory is sufficient at room tem- 
perature, and build a model for bulk crystalline iron 
based upon an spd tight-binding Hamiltonian and 
Stoner theory, treating only the mean atomic spins. 

We therefore identify a minimal set of active phys- 
ical degrees of freedom below Tf to consist of the 
net mean spin moment on each atom and single par- 
ticle states constructed from atomic-like d— orbitals. 
To pass beyond perfect crystalline material, we intro- 
duce a local atomic spin-density extension to Stoner 
theory. 

Construction of Energy Functional — Our selected 
degrees of freedom are (1) the linear combination co- 
efficients 'ipnkih m) describing the bonding among the 
ruz = —2, . . . , 2 atomic d— states of each atom i for 
each band n and each point in the Brillioun zone fc, 
(2) the net spin n„{i)^ a — ±1, associated with each 
atom and (3) the location of each atom i. 

To describe the bonding contribution to the energy 
e„fc associated with each single-particle orbital, we 
use an orthogonal two-center tight-binding Hamilto- 
nian description [p^ . In the usual fashion, we take the 
diagonal elements of the matrix of hopping integrals 
among eigenstates for two atoms separated by a 
distance r along the z-axis, ddb^ dd?:, dda, ddn, ddS, to 
decay exponentially with distance, dd\ — ddXoC^''^ ■ 
To set the primitive matrix elements dd\o and the 
decay factor q we insist that the tight-binding Hamil- 
tonian reproduce well the ah initio density of states 
for the d— bands, which is accomplished by setting 
the Slater-Koster parameters to ddao '■ ddi^o '■ dddo=- 
1 Ryd:l Ryd:-0.5 Ryd and the decay factor to q — 
0.63872. 

In summing over the single particle energies, we 
replace the integral over the Brillioun zone with a 
discrete sum with appropriate weights Wfcp^. The 
tendency of our system to spin-polarize requires us to 
consider separate Fermi occupation numbers fnka for 
each spin channel a. We also include the Fermionic 
entropy 

Snka = -ks [fnka In fnka + (1 - fnka) ln(l - fnka)] , 

so that our band structure energy appears as the first 
sum in the energy functional (|l|) where the enk are the 
eigenvalues of the tight-binding Hamiltonian. All of 
the results below are computed at T = 300°K. 

Combining the Fermi occupations fnka with the 
expansion coefficients ipnka yields the local atomic 



spin densities, n^{i) = Y.n,k,m'^kfnka \i^nk{i,m)\ . 
The Stoner theory of itinerant ferromagnetism de- 
scribes an energetic benefit of polarization due to ex- 
change of magnitude where Nat is the to- 
tal number of atoms in the crystal, m is the spin po- 
larization per atom and / is the Stoner parameter [p^. 
To extend this theory beyond homogeneous bulk sys- 
tems, we introduce a local approximation to the ex- 
change energy in the same spirit as the local-density 
approximation of density functional theory Q. 
In particular, for an inhomogeneous system, we asso- 
ciate a separate energy contribution to the exchange 
from each atom equal to what we would expect on a 
per atom basis from a homogeneous system consisting 
of atoms in the identical environment with identical 
spin polarization. This contribution appears as the 
second sum in our energy functional (|l|) . The physical 
motivation for this approximation is that the itiner- 
ant nature of magnetism in iron tends to smooth vari- 
ations in the spin polarization, limiting the effects of 
gradient corrections. This approach has the advan- 
tage of allowing us to draw upon ab initio values of 
the Stoner exchange parameter in bulk. 

Krasko[^ has performed ab initio linear response 
theory calculations of the Stoner parameter / in bcc 
and fee bulk iron and found it to have a mild, approxi- 
mately linear volume dependence, dl/ds = —0.01 Ry- 
dberg/Bohr where s is the Wigner-Seitz (WS) ra- 
dius, and to have slightly different values for the bcc 
and fee lattices, 1°^^ = 0.072 Ryd//i| and 1°^^ = 
0.069 Ryd/fi% at s = 2.66 Bohr. The use of these 
values in our model gives the correct magnetic and 
non-magnetic states for the bcc and fee structures, 
respectively. However, to yield the correct total en- 
ergy ordering of states, we have found necessary a 
slight enhancement of the bulk Stoner parameters 
to 1°^^ = 0.077 Ryd/At| and If^^ = 0.070 Ryd//i|, 
which leads to no magnetization in the fee phase and 
a bcc phase magnetic moment of 2.56/XB/atom, some- 
what enhanced relative to the accepted moment of 
2.2iiB/atom. 

Going beyond bulk to inhomogeneous systems with 
intermediate coordinations Q, we make a linear in- 
terpolation for li between what would be expected 
at the same WS radius for the Qtcc — 8 and the 
Qfcc = 12 lattices. We determine the local coordi- 
nation number Qj and WS radius Si for each atom 
i with the formulation developed by Sawada|^^ and 
the conversion Si = Ri{a + h/Qi + c/Q'j)/2 from his 
parameter Ri. We find that setting Sawada's coeffi- 
cients to Ai = 4.5023 Bohr-i, A2 = 10.6376 Bohr'^ 
and using a = 1.7144, b = -9.0948, c = 56.372 re- 
produces to within 0.01% the correct coordination 
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numbers and to within 0.3% the correct WS radii for 
the diamond structure, bcc and fee lattices packed at 
the atomic density of bcc iron. With the Stoner pa- 
rameters thus determined, we set the filhng Nd of the 
manifold of d— states, so that the Fermi level for spin- 
down electrons in the bcc structure falls precisely at 
the minimum of the pseudo-gap in the tight-binding 
density of states, reproducing the physical behavior 
observed in ab initio calculations. The resulting fill- 
ing, Nd — 6.7 electrons/atom, is in good agreement 
with the value of Nd — 7.0 used successfully in p^ . 

Finally, for the energy associated with the locations 
of the atoms, we take a power-law relationship be- 
tween the interatomic potential and the hopping ele- 
ments, a standard successful practice in tight-binding 
calculations!^, to produce the final term in (|l|). We 
fit the two parameters b and p to the experimental 
equilibrium lattice constant and bulk modulus for 
the bcc phase of iron, yielding p — 2.2355 Bohr~^ 
and b = 872.5174 Ryd, respectively. The ratio 



X = p/q corresponds to a normalized hardness 21 1, 
a/i = (A — 1)/A « 0.7, in line with the values near 
two-thirds observed previously in tight-binding de- 
scriptions of the transition metals |2|] . 
Our final energy functional is thus, 

= mm' < Wk {fnka^nk - TSnka) 

4' J ^ 

\ nka 
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The constraints on the minimization are Fermi statis- 
tics, < jnka < 1 and the total number of 
d— electrons, = J2nk<7 '^'^f^ka- This formulation 
is equivalent in bulk systems to the familiar formu- 
lation of the Stoner theory in terms of a rigid shift 
between the up and down electronic density of states. 
The present formulation, however, has the advantage 
in treating complex structures of allowing distinct lo- 
cal Stoner parameters to be applied to each atom 
according to its environment. Finally, stationarity 
of the energy functional with respect to the fillings 
dramatically simplifies the evaluation of forces. 

Verification — To confirm the applicability of our 
description to iron, we discuss briefly the comparison 
of our results with available ab initio and experimen- 
tal information for the bulk crystalline phases and for 
grain boundaries in bcc iron. 

In bulk, we reproduce the correct sequence bcc- 
hcp-fcc of phases finding Ehcp—Ebcc=2.2 mRyd/atom 



and Efcc — Ebcc—G-5 mRyd/atom. Our fcc-bcc en- 
ergy difference, is on the order of what is found in 
other calculations |jl^, and theoretical and experi- 
mental extrapolation p3| , p4| . In agreement with 
ab initio calculations pq, |26[ , we observe that the 
hep phase is more stable than the non-magnetic fee 
phase for all values of the WS radius. We predict a 
pressure-induced phase transition from the bcc to the 
hep phase at a WS radius of s = 2.6 Bohr, in good 
agreement with the ab initio studies of |24, 2^. (Note 
that we did not fit our parameters to produce the pre- 
vious two properties.) The small energy differences 
among these phases opens the question of mechani- 
cal stability. Our Hamiltonian gives a mechanically 
stable ferromagnetic bcc phase, even along the Bain 
transformation (C > 0). Our predicted C" and C44, 
which we have made no attempt to fit, are about 25% 
lower than observed experimentally, corresponding to 
an underestimation of about 12% in the frequencies in 
the long- wavelength portion of the phonon spectrum. 

The literature presents an experimental determina- 
tion of an average typical grain boundary energy in 
a-iron and ab initio results for the spin moment dis- 
tribution of the S = 5(310) and E = 3(111) bound- 
aries and the atomic relaxation of the S = 3(111) 
boundary in iron. Table I summarizes our results for 
these and two other symmetric tilt boundaries. Our 
calculations were carried out in supercells contain- 
ing two oppositely oriented boundaries separated by 
at least eighteen layers of atoms. We performed full 
structural and supercell relaxations of these bound- 
aries. 

Table ^ shows that our energy results are in good 
agreement with the experimental studies, which set 
the mean typical grain boundary energy to be ap- 
proximately 770 erg/cm^g. The S = 3(112) bound- 
ary, the coherent twin, is a known special case which 
is expected to be unrepresentatively low in energy. 
The magnetic and structural predictions of our en- 
ergy functional are in excellent quantitative agree- 
ment with the ab initio predictions for the outward 
structural relaxation of the atomic planes of the 
S — 3(111) boundary and for the fractional enhance- 
ment of the spin moments of the symmetry plane for 
both the E = 5(310) and E = 3(111) boundaries. 
(See Az and Am in Table |l|). 

Finally, figure |l| presents a more detailed compar- 
ison with ab initio calculations, showing the spatial 
distribution of spin moments in the vicinity of the 
E — 5(310) boundary. The figure shows that our 
model not only reproduces the enhancement of mo- 
ments on the symmetry plane but also predicts cor- 
rectly the tendency for the spin to fall below the bulk 
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Figure 1: Prediction of atomic spin moments in the 
vicinity of the S — 5(310) symmetric tilt boundary 
in iron. Left: present model, Right: ah initio results. 
Horizontal lines give the bulk phase moment of the 
respective calculation. 



moment before eventually healing back to the bulk 
value as one moves away from the boundary. Some 
discrepancies become apparent in the results at this 
high level level of detail. This comparison serves 
to underscores the fact that there are limitations to 
any simplified, semicmpirical model. Nonetheless, 
the overall level of agreement which we have found 
supports the fundamental soundness of our approach 
and its ability to give accurate predictions of global 
quantities and physical trends for complex structures 
in iron. 

Exchange stabilization of grain boundaries in iron 
— The fact that the embedded atom model consis- 
tently exaggerates boundary energies in iron by a 
factor of two over our functional points to the par- 
ticipation in boundary energetics of a microscopic 
mechanism more subtle than simple distortions of the 
metallic bonding network. The data in Table [| show 
that the exchange interaction is a major factor in the 
physics of the grain boundaries. The tendency to 
minimize the exchange contribution drives the sys- 
tem to lower its energy at the expense of an almost 
compensating dramatic increase in the atomic and 
band structure contributions. The end result of this 
balance is the lower and much more realistic set of 
boundary energies in Table 0. 

Breakdown of the large negative contributions from 
the exchange interaction shows that the enhancement 
comes mostly from an increase by 10-20% of the 
atomic spin moments with a lesser contribution from 



an increase by 2-4% of the Stoner exchange parame- 
ters themselves. (Table |l|.) The increase in magnetic 
moments along the boundary therefore plays the key 
role in the large stabilizing exchange interaction. 

The origin of this increase is the fact that the states 
which contribute to the spin moments in the sys- 
tem tend to be more localized on the grain boundary. 
States with energy below both Fermi levels /i| and m 
are filled equally with up and down spins and so con- 
tribute nothing to the net moment, while those states 
with energy above both Fermi levels are completely 
empty. Only those states between these two Fermi 
levels are filled with uncompensated spins and con- 
tribute to the spin moments. The range of energies 
between the two Fermi levels is near the center of the 
d— band. The states in this energy range tend to lo- 
calize on the boundary as a direct consequence of the 
more open structure of the boundary, which lowers 
the tight-binding matrix elements dd\ and thereby 
narrows the band toward the band-center. We con- 
firmed this latter behavior by direct inspection of the 
electronic states. 

Finally, our model also sheds light on the geomet- 
ric relaxation of the grain boundaries. In all bound- 
aries in our study, the two planes of atoms immedi- 
ately neighboring the symmetry plane to relax out- 
ward and compress into the surrounding bulk. Be- 
cause the Stoner parameters increase with decreas- 
ing atomic volume {dl/ds < 0, above) this dimin- 
ishes the Stoner parameter on the symmetry plane 
but enhances the exchange parameter in a total of 
four planes, one pair on either side of the boundary. 
This relaxation pattern enhances the exchange sta- 
bilization by sacrificing the exchange parameters on 
the single symmetry plane in favor of a total of four 
nearby planes. 

In conclusion, we have developed a local spin den- 
sity functional description of itinerant ferromagnetic 
materials which provides a simple and accurate pic- 
ture of the relationships among geometry, electronic 
structure and stability of tilt grain boundaries in iron. 
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Table 1: Summary of Grain Boundary Results: 
Boundary formation energy [Egb]^ exchange contri- 
bution to the energy [i?ex], change in magnetic mo- 
ment on the boundary plane [Am] , maximum Stoner 
parameter [A/maj;], outward motion of planes imme- 
diately neighboring the boundary [Az]. Results of 
other studies appear in parenthesis. 
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